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Abstract 

We investigated the feasibility of quantum-correlation measurements in nuclear 
physics experiments. In a first approach, we measured spin correlations of singlet- 
spin ( 1 5o) proton pairs, which were generated in 1 H(d, 2 He) and 12 C(d, 2 He) nuclear 
charge-exchange reactions. The experiment was optimized for a clean preparation 
of the 2 He singlet state and offered a 2tt detection geometry for both protons in the 
exit channel. Our results confirm the effectiveness of the setup for theses studies, 
despite limitations of a small data sample recorded during the feasibility studies. 
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1 Introduction 

Entanglement is believed to be a genuine resource for quantum computers and 
quantum communication technology. Entanglement shows up in composite 
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quantum systems where the subsystems do not have pure states of their own. 
This is a strict quantum phenomenon with no classical analogue. Entangled 
states of joint systems are non-local, meaning that the outcome of measure- 
ments performed separately on each subsystem at space-like separation cannot 
be reproduced by local-hidden-variables (LHV) models. Such non-locality can 
be revealed by a violation of an inequality which any LHV model must sat- 
isfy. Such inequality is the Bell-type inequality [1]. Experimental tests of the 
Bell-type inequality have so far been limited to measurements with photons 
[2] rather than measurements with massive Fermions with only one exception: 
a proton-spin correlation measurement performed by Lamehi-Rachti and Mit- 
tig (LRM) about 30 years ago [3]. Note, that quantum non-contextuality has 
recently been tested with massive Fermions in single-neutron interferometry 
experiment by Hasegawa et al. [4]. However, it is well known that, if a theory 
is contextual, it is not necessarily non-local. 



The advantage of using massive Fermions to test Bell-type inequalities is 
that the particles are well localized and the singlet state of the pair can be 
well defined by measuring the internal energy of the two-proton system. In 
this paper, we studied the feasibility of examining spin-correlation measure- 
ments of proton pairs in a 1 So intermediate state generated in 1 H(rf, 2 He)n 
and 12 C(<i, 2 He) 12 B nuclear charge-exchange reactions. By selecting events on 
basis of the structure of the excited state in the remaining nucleus and on 
basis of the internal energy of the 2 He system, we achieve a clean preparation 
of proton pairs under controlled conditions. Our analysis of the experimental 
results described below is compatible with the pioneering LRM experiment. 
However, our experimental setup, the experimental procedure and the data 
analysis improved significantly in comparison to the LRM experiment with 
respect to the following issues: 



a) control of higher order multipole contamination of the singlet state, 

b) control of the contamination due to randomly correlated pairs, 

c) causal separation of the proton pairs, 

d) no preferred quantization axis because of 2ir detection geometry and 

e) record of complete event topology. 



We will structure this paper as follows. In the following section we give a 
brief overview of the experimental arrangement and the data analysis, whereas 
special requirements and the feasibility of spin-correlation measurements are 
emphasized. For details of the experimental setup and the 2 He analysis we 
refer to Ref . [5] . In section 3 we discuss details of the spin-correlation analysis 
and in section 4 our results are presented. 
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Fig. 1. Schematic diagram of the experimental arrangement. Top part: The BBS 
spectrometer in conjunction with the ESN detector system. Lower part: Blown up 
view of the detector setup. The numbers indicate dimensions in units of mm. The 
focal-plane detection system (FPDS) is equipped with two VDCs to determine the 
four-momentum vectors of protons passing the spectrometer. Further downstream, 
the focal-plane polarimeter (FPP) set up, which consists of four MWPCs (Dl, D2, 
D3, D4) and two plastic-scintillator arrrays (SI and S2). A carbon analyzer (C) is 
placed just next to the scintillator SI. 

2 Experimental setup 

The measurements were carried out using 172 MeV deuteron beams provided 
by the AGOR cyclotron of the Kernfysisch Versneller Instituut (KVI), Gronin- 
gen. The deuterons were incident on a carbon foil of thickness 9.4 mg/cm 2 , 
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which was mounted in the scattering chamber of the Big-Bite Spectrometer 
(BBS) [6]. A schematic view of the experimental setup is shown in Fig. 1. 
The EuroSuperNova (ESN) detector system used consists of a pair of gas- 
filled Vertical-Drift Chambers (VDCs) for momentum reconstruction of the 
protons and a focal-plane polarimeter, which comprises Multi-Wire Propor- 
tional Chambers (MWPCs) and a pair of scintillator paddles SI and S2 for 
time-of-flight and energy- loss measurements [5]. 

The 1 S'o proton pairs were prepared in a 12 C(d, 2 He) 12 B nuclear charge-exchange 
reaction. A reaction of type (d,2p) is referred to as (d, 2 He) if the outgoing 
protons couple to the 1 Sq state. The 2 He system is unbound by an internal 
energy of about 0.5 MeV, which is defined by the maximum of the (pp) 1 S'o 
final-state-interaction strength. The proton pairs emerging from the charge- 
exchange reaction were momentum analyzed in the BBS spectrometer, which 
was positioned at an angle Obbs = 0°. The extreme forward angle was chosen 
to minimize the angular momentum transfer in the reaction which favors pure 
spin-flip (Gamow- Teller) type transitions and puts an additional constraint 
on the 1 S'o character of the proton pairs. 



2.1 2 He identification 

Due to the finite momentum and angular acceptances of the BBS focal plane, 
only proton pairs of relative kinetic energies less than 1 MeV in the 2 He center 
of mass were detected as depicted in Fig. 2. The spectrometer therefore acted 
as a highly exclusive filter and contaminations of higher-order multipoles were 
limited to the percent level [7]. A dominant background of randomly correlated 
protons was due to the breakup of the deuteron, a reaction yielding protons 
with a momentum overlapping with the momentum range of interest. A clean 
identification of 2 He events therefore necessitated a proper reconstruction of 
the excitation energy of the residual nucleus 12 B and determination of the 
relative timing of the two correlated protons with a good resolution. 

An event-trigger condition required that at least one proton passed through 
the spectrometer and was registered by coincident signals in scintillator planes 
SI and S2. The coincidence window was set to less than 20 ns in order to 
minimize randoms caused by particles originating from different beam bursts 
each separated by 23 ns. The data, read out from the MWPCs, were fed 
into a fast online processing system, where they were tested on double-track 
conditions. Those events passing the test were stored for offline analysis. 

The momentum vector and the focal-plane interception time of protons were 
determined in offline analysis [5] from the VDC data. After being triggered, the 
VDC TDC channels remained active for about 380 ns, which is the maximum 
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Fig. 2. Typical relative kinetic-energy distribution in the 2 He center of mass obtained 
from the measured momenta of coincident protons. 



drift-time associated with the events. A typical spectrum of the difference 
in focal-plane interception time corresponding to kinetic energies below the 
12 B Q-value is shown in the top left of Fig. 3. The 2 He protons, which can 
for kinematical reasons be separated by at most a few 1CT 9 s (ns) in the 
BBS focal-plane, appear as a dominant prompt peak. The interception-time 
spectrum shows a peak centered at t=0 and satellite peaks due to random 
coincidences. The satellite peaks are due to randomly correlated proton pairs 
reflecting the beam-burst repetition rate. An effective identification of 2 He 
events could be achieved by requiring an interception-time difference in the 
window ±10 ns. An energy spectrum accumulated under this condition is 
shown in the lower right of Fig. 3. The kinetic energies were calculated for 
two-proton events assuming 2 He kinematics [5]. The energy spectrum shows 
two prominent peaks superimposed on a continuous background (see top right 
of Fig. 3). The peak at 169 MeV corresponds to the Q- value of the 1 H(d, 2 He)n 
reaction and is due to a hydrogen contamination of the 12 C target. The peak at 
157 MeV corresponds to the Q-value (Q = -14.81 MeV) of the 12 C(d, 2 He) 12 B 
reaction. The structures at lower kinetic energies correspond to transitions 
to 12 B excited states. Unphysical energies larger than in the incident beam 
energy of 172 MeV are due to randomly correlated protons originating from the 
deuteron breakup. For the spin-correlation analysis, it was further required, 
that the sum of the kinetic energies of the proton pairs was equal to or less 
than the 12 B threshold or in an energy window defined by the position and 
the width of the neutron peak. 
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Fig. 3. Top left: Focal-plane interception-time difference for detected proton pairs. 
For the data set shown the random contamination is about 15%. Top right: Raw- 
kinetic energy spectrum of detected proton pairs. Prominent peaks indicate 2 He 
protons originating from the 12 B ground-state transition and the p-n transition. 
Lower left: Kinetic energy spectrum gated on satellite peaks in the interception-time 
difference (no correlated protons). Lower right: Kinetic energy spectrum gated on 
the prompt peak in the interception-time difference. 



3 Spin-correlation measurements and analysis 



In order to measure the scattering angle in the carbon analyzer, identification 
of proton tracks upstream and downstream of the analyzer was required. We 
followed the fate of each proton [8] as it passed through the carbon analyzer, 
acquiring the information in the detector systems D1-D4 and SI, S2. We used 
only those events where both protons scattered into an angular range larger 
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Fig. 4. Geometry applied for extracting the experimental correlation function Eq. 
3. The primary target is indicated by 12 C, the primary scattering normal n2jj e 
refers to the incoming deuteron and 2 He center-of-mass movement, y is the normal 
to the BBS (horizontal) symmetry plane. The momenta of proton 1 and proton 2 
upstream the analyzer are indicated bypi( 2 ), downstream the analyzer by p^ 2 )> an d 
the normal and angle of the analyzing scattering for proton 1 (2) are indicated by 
n X ( 2 ) and #1(2), respectively. The angle between the proton momenta Vi{2) 15 given 
by 9 pp . 

than 3 degrees in the carbon analyzer, since the most forward scattering is 
predominantly Coulomb type which is not spin dependent. The experimental 
setup provided the flexibility to arbitrarily choose the reference axis during 
the off-line analysis. 

The correlation function P(0) for two spin-| states can be measured according 
to 

= n ++ + m_ - jy + _ - iv_ + = n ++ + n__ - ivy - iv_ + 

U N ++ + N__ + N + _ + iV_ + N total ' 1 } 

where 6 is the angle between two arbitrary quantization directions Qi and 

— * 

Q 2 orthogonal to the momenta pi and p 2 of the two correlated protons, i.e. 

Qi _L pi and Q 2 -L P2 and N ++ (N ) is the number of events, where both 

protons scatter to the left (right) of the quantization direction and iV + _ (iV_ + ) 
is the number of events, where proton 1 scatters to the left (right) and proton 
2 scatters to the right (left) and N tota i is the total number of events. Taking 
into account the finite analyzing power of the carbon analyzer, the number of 
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events, in the above correlation function have to be weighted on an event-to- 
event basis by the analyzing powers Ay> 2 (0 lj2 , E 1>2 ) [9] 



yielding the experimental correlation function P exp (see also LRM 

ATW _|_ T\TW _ ATW _ T\TW 

Pe*M = ^^"^ ^ - + • (3) 

^ total 



The geometry applied to extract the experimental correlation function Eq. 3 
is shown in Fig. 4 and Fig. 5. The correlation analysis was only applied for 
events where the direction of the primary scattering normal Hi deviated less 
than 1° from the normal of the BBS symmetry plane y. Due to this selection, 
our analysis became compatible with the LRM analysis [3]. 

The sign convention for the correlations is shown in Fig. 5, where the example 

of a right-right ( ) correlation is depicted. The quantization axis Qi( 2 ) is 

chosen along the scattering normal nu 2 ) with the convention, that the projec- 
tion on the normal the BBS symmetry plane Qi(2) • V is positive. According to 
this convention, a positive definite projection of the scattering normal n^) • y 
indicates scattering to the left, a negative definite projection scattering to the 
right of the quantization axis. 

The correlation angle 9 (see Fig. 5) is defined as 

e = \(px- <p 2 \ , (4) 

whereas the angle <fi( 2 ) is measured in respect to the vector xu 2 ^ given by 

^1,2 = V X Pl(2) • (5) 



The definition of ipu 2 ) and 6 holds also for finite 6i because of the purely 
transverse character of the analyzing reaction. 

In quantum theory, the operator that corresponds to the correlation function 
is 

P = Qi-a®Q 2 -a , (6) 



acting in the Hilbert space Hi <8> H 2 in 2 ® 2 dimension and a are the Pauli 
matrices. The correlation function Pqm is given by the mean value of this 
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operator. For a pure state this correlation function could be easily computed. 
For a singlet state we have 



P%T(9) = - cos 9. 



(7) 



However, if the state is mixed the mean value should be averaged over the 
ensemble. Taking into account the effect of a random contamination of the 
pure singlet state, the quantum expectation deviates from Eq. 7. In fact, we 
introduce a factor 7 which interpolates between the unpolarized state 1/4 and 
the singlet state = (| + -) - | - +))/yf(2), with / the unit matrix. 
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Fig. 5. Sign convention for the correlations depicted for a right-right ( ) event 



and 6 PP = 0°. 



The density matrix of such a state is called Werner states [10] and it is given 
by 



p w = (1-7)4+7^ )<*" 



(8) 



This effect reduces the quantum expectation value of the correlation functions 
as follows 



pWerner {9) = _ lcos9 



(9) 



3 
2 
1 



-1 
-2 
-3 

20 40 60 80 100 120 140 160 180 
8[deg] 

Fig. 6. The extracted experimental correlation function P exp {9) (Eq. 3) in compar- 
ison with the quantum mechanics prediction P^Q M rner (Eq. 9) for 7 = 0.9 

4 Results 

During two days of data taking, the mean rate of 2 He events ending up in the 
correlation analysis amounted to about 0.1 Hz/nA, despite the fact that the 
(d, 2 He) production rate was nearly a factor 50 higher. This loss in statistics was 
mainly due to inefficiencies in the particle tracking downstream the analyzer. 
For the final analysis the data were binned into 12 angular 9 bins, which 
yielded on average 10 3 events per angle bin. 

The random events, i.e. events shown in the lower left of Fig. 3, yielded an 
averaged correlation P exp = 0.05 ±0.7 per angle, which justifies the treatment 
of the background as a random contribution as indicated in Eq. 9. The overall 
mean of random events contributing to the prompt events amounted to 10%. 

In Fig. 6 we show the extracted experimental correlation function P exp in 
comparison with the quantum mechanics prediction for Werner states Eq. 9 
for 7 = 0.9. Despite of the large uncertainties the data exhibit a trend which 
agrees well with the quantum-mechanics prediction and yields a x^/d.o.f = 
Yji{Plx P ~ Pqm) 1 1 ^Plxp) 2 = 0-96, which has to be compared with x 2 /d.o.f = 
1.93 if P QM is replaced by P const = 0. 
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Fig. 7. CHSH inequality for all possible combinations found from the probability 
shown in Fig. 6. The dashed line shows the experimental results, the dotted line the 
QM prediction. The shaded band represents the statistical error. 



In order to demonstrate the power of the presented experimental approach, 
we further used the values we obtained for P exp to extract the values for a 
correlation function proposed by Clauser, Home, Shimony and Holt (CHSH) 
[11] which can be written in the form 

\P(0 a , 9 b , 9 C )\ = \P{6 a ) - P(9 a + 9 b )\ + \P(0 a + 9 b + 9 C ) + P(9 a + 6 C )\ , (10) 



with # a ,b,c various sets of correlation angle 9. All possible angular combinations 
yield 84 test cases which have been plotted in Fig. 7 in comparison with the 
quantum-mechanics predictions. A discussion of the CHSH-type correlations 
is beyond the scope of the present paper and will be the topic of a forthcoming 
publication [12]. From Fig. 7 it becomes obvious that the present experimental 
data set suffers from large statistical uncertainties. Nevertheless, the data 
exhibit a tendency to stay below the quantum mechanical results. A fact 
which would point towards classical scenarios. 
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5 Conclusion 



In this paper, we presented a new experimental approach to study the feasi- 
bility of examining spin-correlations measurements in nuclear physics. With 
an improved detector setup, which removes the ambiguity in the track recon- 
struction, measurements with significant precision will become feasible. We 
are convinced that our experiments will have many potential applications in 
future quantum communication technology. 

This work was performed as part of the research program of the Stichting voor 
Fundamenteel Onderzoek der Materie (FOM) with financial support from the 
Nederlandse Organisatie voor Wetenschappelijk Onderzoek . It was supported 
by the NSERC Canada, the European Union through the Human Capital and 
Mobility Program and the Fund for Scientific Research (FSR) Flanders. 
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